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Abstract 

All the possible super-conducting order parameters for the LaOFeAs system are classified by their 
transformation under the complete crystal symmetry. The general forms of the super-conducting 
gap functions for each class are discussed. We find that the gap functions in such a multi-band 
system belong to three types, full gap, nodal type and finite "Fermi arc" type. Possible physical 
consequences caused by different types of gap functions are also discussed. 
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I. INTRODUCTION 



The iron based superconductors, the second family of the high temperature supercon- 
ducting materials, have been discovered very recently. Hosono group first obtained super- 
conductivity in LaOFeAs with = 2QK by replacing some Oxygen atom with F jl|. Soon 
after it, super-conductivity is discovered in several similar materials with Tc = 41 K in 
CeOi-xFxFeAs [2|] and Tc = 43 K in SmOi-xFxFeAs ^. Further replacing La atoms with 
other rare earth elements rapidly raises Tc up to b2K in Pr[Oi_xFx]F eAs {4] and 557^ in 
Sm[0 i-xF^FeAs [5||. With some other ways to introduce carriers, super-conductivitv with 
the T(.= 55K in Gdi^xThxO F eAs and ReFeAsOi^s are obtained by several groups p, 5, ^■ 

N. L. Wang's and Z. Fang's group first reported that the ground state of the parent 
compound shows SDW long range order with stripe like spin configuration by combining 
the first principle calculation and optical measurement [9|. The SDW state has been con- 



firmed by neutron scattering results from two independent groups 
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calculations 



tion at 135K M 
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ll| . and several other 



13l|. The NMR data also show strong evidence of magnetic phase transi- 
15l |. which is consistent with the previous neutron scattering and optical 



conductivity experiments. Upon F doping, the SDW is suppressed very quickly and super- 
conductivity appears. Since the maximum Tc achieved in this family is around 55K, which 
is well above the McMillan limit, the super-conductivity here is likely to be non-BCS type. 
The pairing symmetry of the super-conductivity in iron based super-conductors is one of 
'he strong mag netic fiuctuation, which is suggested by the first principle 



the key issues, 
calculation 
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strong mag net] 



and NMR measurement 
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15| . is considered by many 



authors to be the "pairing glue". However it is still under debate if the ferromagnetic 
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22| or anti- ferromagnetic [12, ll3| fluctuation is responsible for the super-conductivity 



here. The ferromagnetic (anti-ferromagnetic) fluctuation will mediate attractive interaction 



in the spin trip 



state 
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et (singlet) channel and thus leads to the triplet (singlet) super-conducting 
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28l |. The current experimental data still can not confirm the spin 



state of the Cooper pairs and band structure calculation indicates nearly degenerate multiple 
Fermi Surfaces (FS) in the system [22||. Therefore at the current stage, there still exist many 
possible pairing order parameters leading to many different types of gap functions . It 
is then important to classify all these possible super-conducting order parameters by their 
different transformation behavior under the full symmetry of the crystals, which will be done 
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in the following sections. After the classification, we will obtain many classes of the super- 
conducting order parameters with each of them forming an irreducible representation of the 
full symmetry the system have. Also we will discuss the general form of the gap functions 
for each class of the super-conducting phase, which can be grouped into three catalogs. The 
simplest type, which is referred to type (I) gap function, has a fully opened gap over the 
whole FS for any non-zero order parameter. The type (II) gap function contains a nodal 
point on the FS where it vanishes. And the type (III) gap function will vanish on a finite 
section of the FS, which is mentioned as "Fermi arc" in the previous literature [23 1. 



II. CLASSIFICATION OF THE SUPER-CONDUCTING ORDER PARAMETERS 

The full symmetry of the crystal can be view as the direct production of time reversal, 
spacial inversion, charge conservation, translational symmetry and the point group symmetry. 
In the present study, we assume that the appearance of super-conductivity only breaks the 
charge conservation and keeps all the other symmetry. The purpose of the current paper is 
to classify all the possible super-conducting order parameters by the group representation 
theory. 



A. The Fermi surface topology and the point group symmetry of LaOFeAs 



In the absence of SDW state (after F-doping), LaOFeAs crystallizes in layered square 
lattice with space group P4/nmm. There are two Fe sites per unit cell, which are coordinated 
by As-tetrahedrons. The point group symmetry of LaOFeAs is ^4^, which contains a 4- 
fold rotational axis along z-direction, four vertical reflection planes along x, y, xy and xy 
directions, a horizontal reflection plane, as well as the spacial inversion. According to the 
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18|, 
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20I, 



2l|there are two 



first principle calculation, for non magnetic LaOFeAs, |16l . 
hole type Fermi surfaces enclosing F point and two electron type Fermi surfaces enclosing M 
point. Upon F doping, the three hole type Fermi surfaces shrink rapidly and will be neglected 
in the present paper. Therefore we can focus on the super-conducting order parameters on 
the two electron type Fermi surfaces. The shape of the Fermi surfaces enclosing M point is 
very close to ellipse with their long axis along (1, 1) or (1, 1) directions. The first principle 
calculation indicates that the states around the FS around M point are mainly from the 
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Figure 1: The schematic Fermi surface of the LaOFeAs, the FS of band 1 (2) is plotted by the black 
(red) ellipse. 



xy,xzand yz iron 3d orbitals with very strong angle dependence along the FS [22|. Before 
we begin to classify all the possible superconducting order parameter A (k)for a multi- 
band system, we first have to define the Bloch states (r) around M point properly 
to guarantee'0^^ (r) to be smoothly connected along the momentum space. This "smooth 
condition" give us unique definition of the energy bands, which gives the ellipse FS for band 
"1" ( "2") along (1, 1) ((1, 1) )directions as shown in Fig. 1. 

By definition the origin for all the point group operations and spacial inversion should be 
located at F point. While one can easily prove that using the periodicity of the momentum 
space, the origin can be moved to M point. Therefore in the following of the paper, all the 
point group and spacial inversion operations are taken respect to the M point. 

With such a smooth condition for the energy bands, we can write down the effect of the 
-D4/1 group elements gi acting on the Bloch functions in the two bands forming the FS around 
the M point. In general it can be written as 



aii'nk ('^) = DZn {91) ^„DR(-)(g,)fc (^) 

where 2x2 matrices D'j^^ [gi] are the representation matrices in the orbital space, which 
equal to unit matrix for operators (i?, C2 , 2cr^)and for the rest of them. The above 
representation in the orbital space can be checked by the symmetry right on the FS. From 
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Fig.l, one can easily find that the FS of band "1" ("2") remains unchanged under the first 
catalog of the operators and transforms to FS of band "2" ("1") under the second catalog of 
the operators. 



B. The super-conducting order parameter and its transformation under D^^ 

In the multi-band system, the super-conducting order parameter can be written as 

where a, P = 1,2 are the band Indies, //, =|, t denote the electron spin. The functions 
A°f^ (k) form a representation of the full symmetry of the system, which can be expressed 
as the direct production of time reversal, spacial inversion and point group symmetry. In 
addition, the anti-symmetric nature of the fermionic many body wave function requires 
A^^ (fc) = — A^^ (~^)- In the present study, we follow the treatment of Sigrist and Ueda 



on the effect of spin-orbital coupling (SOC)[29||, assuming the SOC is not too strong so that 
the induced coupling between the singlet and triplet pairing states can be safely neglected. 
Therefore we can first divide A°f^ (k) into singlet and triplet two subgroups. Among them, 
the singlet pairing order parameter can be written as a 2 x 2 matrix A'^ingiet (^) ^'^^^ the 
following requirement from the anti-symmetric condition upon particle exchange, 

^singlet (^) ~ '^singlet (^^) 

. Therefore for even (odd) parity state A"^^ {k) must be a symmetric (anti-symmetric) 



matrix. Following the standard notation the order parameters for the triplet pairing statefSOfl 
can be expressed by a vector '^"'^ (A;) with each component of the d- vector to be a 2x2 matrix. 
And the requirement from the anti-symmetric condition here is 

, where for even (odd) parity states '^"'^ (A;) must be anti-symmetric (symmetric) matri- 
ces. With the above notation of the super-conducting order parameters for both singlet 
and triplet pairing states, we can write down the symmetry transformations of these order 
parameters under time reversal and point group operators. For time reversal symmetry we 
have 
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, for spin singlet states and 

At last we consider the symmetry transformation under point group operators, which are 
group elements of D4/J. As we have mentioned before, in the present study, the energy scale 
of spin-orbital coupling is considered to be small compared with the energy splitting between 
the spin spin singlet and triplet pairing states but is still big enough to lift the degeneracy 
in the spin space spanned by vector As the consequence, the spin degree of freedom 

of the Cooper pairs has to be frozen to the spacial motion of the electrons. Therefore the 
d-vector has to rotate exactly following the momentum k under the point group operators, 
which leads to 

gi'^^P (k) = {gi) (^«(-) {gi) k) 

where l)^(='=) is the representation in three-dimensional space with positive (for spin space) 
or negative (for k-space) spacial inversion operation respectively. 

By the transformation under symmetry operators, we can construct a representation of 
the full crystal symmetry using the super-conducting order parameters as the basis. As we 
have mentioned above, the two corresponding Bloch bands we are interested in are defined 
smoothly in k-space. Therefore the SC order parameters A^^ [k) should be also smooth 
function in k-space, which allows us to expand them by first three sphere harmonics and use 
them as the basis for k dependence. 

In general such a representation is reducible and can be reduced to a serial of irreducible 
representations. In the following two sections, we will list all the irreducible representa- 
tions and discuss the corresponding basis functions for both the singlet and triplet super- 
conducting states in LaOFeAs. 
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Table I: Basis gap functions for the symmetry D^h (a). Spin-Singlet 
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C. The irreducible representation and basis gap functions for spin singlet super- 
conducting state in LaOFeAs. 

In the above table, we list all the basis gap functions for spin-singlet super-conducting 
states in LaOFeAs, where (Jo , en, cr2 , C3 are 2x2 unit matrix and Pauli matrices respec- 
tively, which form a complete basis for the orbital space. On the first three columns we list 
the even-odd properties of the basis function under the time reversal (T), spacial inversion 
(P), spin exchange {S) and orbital exchange (O), with the plus and minus sign represent- 
ing the even and odd functions under the transformation respectively. The anti-symmetric 
condition for the Fermion exchange requires P ■ S ■ O — —1, which is well satisfied by each 
basis function. Due to the multi-band nature, the system show some interesting feature 
which is different with the single band case. Firstly because of the additional orbital degree 
of freedom, it is possible to have odd parity basis functions for spin singlet phase, as in 
representation and E^. Secondly in Aig and i?2g, the s-wave basis functions can mix 
with some d-wave basis functions. This is because the special FS topology around the M 
point. The two electron type FS enclosing M point form a 2-dimensional representation for 
Dih. If we only consider the pairing order parameters within one particular band, the D^h 
symmetry will be broken. 

In the following, we will discuss in deal the super-conducting gap functions for Aig, which 
is the trivial representation of D/^h, as an example. In this case, the super-conducting order 
parameter can be written as the linear combination of the three in-plane basis functions. 



7 



which reads, 




(A;) = a {cos{kx) + cos{ky)) + hsin{kx)sin{ky) 



/S.g (k) = a{cos{kx) + cos{ky)) — hsin{kx)sin{ky) 




{k) = c {cos{kx) + cos{ky)) 



With the different value of a, 6 and c, the super-conducting gap may have following three 
different types. 

1. Full energy gap, which is labeled by I in the above table. When the intra-band s-wave 
component a 7^ 0, there is finite energy gap right on the Fermi level for the whole FS. 

2. Energy gap with a nodal line along the x and y axises, which is labeled by // . In 
order to have nodal lines in this particular representation, we must have a = c = 
and 6 7^ 0. In this case the super-conducting gap vanishes at four FS crossing points. 

3. Fermi arc which is labeled by ///. With this type of gap structure, there is a finite 
section of FS which remains ungaped, which is first proposed by us in reference [23 1 . 
This situation may happen when only inter-band pairing strength c 7^ 0. 

The second important representation for the singlet pairing states in LaOFeAs is , 
within which the SC order parameters are linear combinations of the intra-band and inter- 
band (ia,2_j^2 states. In this case the pairing is the strongest at four FS crossing points and 
vanishes along the x and y axises. When only the inter-band component is non-zero, the 
gap function will have a "Fermi arc". 

D. The irreducible representation and basis gap functions for spin triplet super- 
conducting state in LaOFeAs. 

In the above table, we list all the basis gap functions for spin-triplet super-conducting 
states in LaOFeAs. Since the total spin of the Cooper pairs generate extra degree of 
freedom, the situation for the spin triplet SC states is more complicated. In the present 
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Table II: Basis gap functions for the symmetry D^h (b). Spin- Triplet 
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study, we will only focus on the possible gap structure. The first five representations with 
even parity are purely inter-band pairing states. In the previous paper, we have discussed 
the gap structure of inter-band s-wave and dx^-yi-wave states , which are listed in theSig 
and Aig representations respectively. For the inter-band s-wave state with order parameter 
^0^2 [62(0;^ + y"^)], the system will have full gap at the FS when the amplitude of the pairing 
order parameter is much bigger than the maximum splitting of the two bands at the FS 



and will have "Fermi Arc" behavior otherwise 



231 . And as we have discussed in the previous 



paper, for the inter-band dx2-y2 state, the "Fermi Arc" will always appear around the diagonal 
line. The basis functions in thei?2g and A2g representations describe the inter-band pairing 
states along the c-axis, which are very unlikely for a layered compound like LaOFeAs. 

The basis functions in the rest of the triplet representations mainly describe the intra- 
band p-wave pairing states. But they can be mixed with the corresponding inter-band p- 
wave pairing states. The Aiu representation contains three in-plane pairing states, namely 
<^o[(^x{x) + ey{y)], (J^lexiy) + ey{x)] and aile^ix) + ey{y)]. The first one is quite well known 
in the Helium III system as the BW phase. While in the present two-band system, very 
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interestingly the BW phase here may be mixed with two other states. The first one is an 
intra-band pairing state with the order parameter ~^ — (J3[ex{y) + ey{x)], and another one is 
an inter-band pairing state with ~^ — ai[ex{x) + ey{y)], which is nothing but the inter-band 
BW phase. Although the intra-band BW phase itself has a full energy gap on the FS, it 
may acquire line nodes or "Fermi Arc" after mixing with the intra-band BW phase. We will 
discuss the physical properties of this state in detail in a separated paper. The behavior of 
the energy gaps for the A2u,Biu and states are quite similar with that of Aiu states. 
Another important representation for the triplet states is E^, which contains the ABM phase 
in Helium III system. The intra-band ABM phase, which is proposed by M. Sigrist and M. 
Rice for the super-conductivity in StRuO^, can be obtained by linear combination of two 
basis, namely a"o[e2(a;)] and ao[ez{y)]. The intra-band ABM phase breaks the time reversal 
symmetry and also can mix with the inter-band ABM phase. Similarly for purely inter-band 
ABM phase, the "Fermi Arc" behavior will appear. 



III. CONCLUSIONS 



In the present study, we classified the possible super-conducting order parameters for 
LaOFeAs systems respect to all the crystal symmetry including time reversal, inversion 
and point group symmetry. Since the spin-orbital coupling in the system is weak, we can 
still divide the possible super-conducting order parameters into spin singlet and triplet two 
classes. The for each class we give the detailed classification table with each class forming 
an irreducible representation. We also discussed the possible gap structure for each class of 
super-conducting phases. In the LaOFeAs system, partly due to its multi-band nature, there 
are three possible types of gap functions, namely full gap, nodal gap and "Fermi Arc" type. 
The different types of super-conducting gap will have very different thermal dynamic and 
super-conducting properties in low temperature, i.e. the specific heat and superfiuid density. 
The super-conducting state with full gap, like s-wave, will have exponential temperature 
dependence in low temperature for specific heat and superfiuid density. If there are line nodes 
on the FS, both the specific heat and superfiuid density will show power law temperature 
dependence. While for super-conducting state with finite "Fermi Arc", the low temperature 
behavior of specific heat appears like there are some "residue density of states" on the 
Fermi level. Currently there are some experimental evidence indicating the existence of low 
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31 



32| . But we still need further 



energy quasi-particles in the super-conducting phase(l5. 
experimental data to make the conclusion. 

We acknowledge valuable discussions with Y. P. Wang, N. L. Wang, J.R. Shi, F.C. Zhang, 
Lu Yu and Z.B. Su. This work is supported by NSF of China and 973 project of the Ministry 
of Science and Technology of China.. 



[1] Y. Kamihara, et. al., J. Am. Chem. Soc, (doi:10.1021/ja800073m), (2008). 

[2] G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, P. Zheng, J. L. Luo, and N. L. Wang, 

cond-mat/0803.3790vl, (2008). 
[3] X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen and D. F. Fang, cond-mat/0803.3603vl, 

(2008). 

[4] Z. A. Ren, J. Yang, W. Lu, W. Yi, G. C. Che, X. L. Dong, L. L. Sun, Z. X. Zhao, cond- 
mat/0803.4283vl, (2008). 

[5] Zhi-An Ren, Wei Lu, Jie Yang, Wei Yi, Xiao-Li Shen, Zheng-Cai Li, Guang-Can Che, Xiao-Li 
Dong, Li-Ling Sun, Fang Zhouand Zhong-Xian Zhao, cond-mat/0804.2053, (2008). 

[6] Cao Wang, Linjun Li, Shun Chi, Zengwei Zhu, Zhi Ren, Yuke Li, Yuetao Wang, Xiao Lin, 
Yongkang Luo, Xiangfan Xu, Guanghan Cao, and Zhu'an Xu, cond-mat/0804.4290, (2008). 

[7] Zhi-An Ren, Guang-Can Che, Xiao-Li Dong, Jie Yang, Wei Lu, Wei Yi, Xiao-Li Shen, Zheng- 
Cai Li, Li-Ling Sun, Fang Zhou, and Zhong-Xian Zhao, cond-mat/0804.2582, (2008). 

[8] H. H. Wen, and et.al., cond-mat/0803.3021vl (2008). 

[9] J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. Wu, G. F. Chen, X. Dai, J. L. Luo, 

Z. Fang, and N. L. Wang, cond-mat/0803.3426vl, (2008). 
[10] Clarina de la Cruz, Q. Huang, J. W. Lynn, Jiying Li, W. Ratcliff H, J. L. Zarestky, H. A. 

Mook, G. F. Chen, J. L. Luo, N. L. Wang,and Pengcheng Dai, cond-mat/0804.0795, (2008). 
[11] M. A. McGuire, A. D. Christianson, A. S. Sefat, R. Jin, E. A. Payzant, B. C. Sales, M. D. 

Lumsden and D. Mandrus,cond-mat/0804.0796, (2008). 
[12] Fengjie Ma, Zhong-Yi Lu and Tao Xiang, cond-mat/0804.3370, (2008). 
[13] T. Yildirim, cond-mat/0804.2252, (2008). 

[14] K. Ahilan, F. L. Ning, T. Imai, A. S. Sefat, R. Jin, M.A. McGuire, B.C. Sales, D. Mandrus, 
cond-mat/0804.4026, (2008). 



11 



[15] Yusuke Nakai, Kenji Ishida, Yoichi Kamihara, Masahiro Hirano, Hideo Hosono, cond- 
mat/0804.4765, (2008). 

[16] G. Xu, W. Ming, Y. Yao, X. Dai, S. C. Zhang, and Z. Fang, cond-mat/0803.1282v2, (2008). 
[17] D. J. Singh and M. H. Du, cond-mat/0803.0429vl, (2008). 

[18] I.I. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du, cond-mat/0803.2740.vl, (2008). 
[19] K. Kuroki, and et.al., cond-mat/0803.3325vl (2008). 
[20] S. Lebegue, Phys. Rev. B 75, 035110 (2007). 

[21] K. Haule, J. H. Shim, and G. Kothar, cond-mat/0803.1279vl (2008). 
[22] H. J. Zhang, et. al. cond-mat/0803.4487. 

[23] X. Dai, Z. Fang, Y. Zhou. F. C. Zhang, cond-mat/0809.3982vl (2008). 

[24] Xiao-Liang Qi, S. Raghu, Chao-Xing Liu, D. J. Scalapino, Shou-Cheng Zhang, cond- 

mat/0804.4332 (2008). 
[25] Patrick A. Lee, Xiao-Gang Wen, cond-mat/0804.1739 (2008). 
[26] Zi-Jian Yao, Jian-Xin Li, Z. D. Wang, cond-mat/0804.4166 (2008). 
[27] Qimiao Si, Ehhu Abrahams , cond-mat/0804.2480 (2008). 
[28] T.A. Maier, D.J. Scalapino, cond-mat/0805.0316 (2008). 
[29] M. Sigrist and K. Ueda, Rev. Mod. Phys. 63, 239 (1991). 
[30] A. J. Leggett, Rev. Mod. Phys., 47, 331 (1975). 

[31] H. Yang, et ah, cond-mat/0803.0623; G. F. Chen, et ah, cond-mat/0803.0128 (2008). 
[32] L. Shan, et ah, cond-mat/0803.2405v2 (2008); X. Zhu, et ah, cond-mat/0803.1288 (2008); G. 
Mu, et ah, cond-mat/ 0803.0928 (2008). 



12 



